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Abstract 

Using newly available form factors obtained from light cone QCD sum rules in 
full theory, we study the flavor changing neutral current transition of — > 
decay in the family non-universal Z' model. In particular, we evaluate the differential 
branching ratio, forward-backward asymmetry as well as some related asymmetry 
parameters and polarizations. We compare the obtained results with the predictions 
of the standard model. It is demonstrated that most of the mentioned observables are 
sensitive to family non-universal Z' gauge boson. The order of differential branching 
ratio shows that this decay mode can be checked at LHC in near future. 
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1 Introduction 



The heavy baryons containing a single heavy quark constitute a perfect laboratory to test 
the non-perturbative aspects of QCD. The A b — > A/i + yU~ transition which is based on the 
flavor changing neutral current (FCNC) transition of b — > sfi + fi~ at quark level has been 
recently observed in CDF collaboration at FermiLab [1]. It is also planned to be checked 
at LHCb collaboration at CERN [2]. 

The theoretical studies on the branching ratio of £& — > Yjt + t~ in standard model (SM) 
[3] show that this decay mode is also possible to be observed at LHC. It is expected that 
the experimental studies on the heavy baryons and their decay properties constitute one 
of the main direction of research program at LHC. Hence, the theoretical calculations can 
play an essential role in this regard. 

Although the SM has predictions in perfect agreement with collider data up to now, 
there are some problems such as neutrino oscillations, baryon asymmetry, unification, dark 
matter, strong CP violation and the hierarchy problem, etc. which can not be addressed 
by the SM and still remain unsolved. To cure these deficiencies, there are a plenty of new 
physics (NP) models such as different extra dimension models (ED), various supersymmetric 
(SUSY) scenarios, etc. One of the most important new physics scenarios is Z' model, 
appears in many grand unified theories, such as SU(5) or string-inspired E6 models [4- 
8]. The two Z' models in agenda are family non-universal Z' [9,10] and leptophobic Z' 
scenarios [11, 12]. 

The idea of extra heavy Z boson comes from the extension of gauge group SU(5), 
predicted by the grand unification theories to larger group SO (10). The 5*0(10) gauge 
group is the next important one after SU(5) having one extra rank. Hence, this gauge 
group requires at least one extra neutral gauge boson [13]. In general, Z' gauge couplings 
are family universal [10,14-19], however, due to different constructions of the different 
families, in string models it is possible to have family non-universal Z' couplings. In some 
of them, three generation of leptons and also the first and second generation of quarks have 
different coupling to Z 1 boson, when compared to the third families of quarks [10,20,21]. 
For more information about this model see for instance [9, 10,22-26]. The leptophobic Z' 
model means that the new neutral gauge boson does not couple to the ordinary SM charged 
leptons. 

The study of the Z' phenomenology is an important part of the scientific program of 
every present and future colliders and due to its heaviness, Z 1 may be used to calibrate 
future detectors [13]. For constraints on the mass of Z' boson and mixing parameters of 
this model see for example [27-29]. There are direct search for Z 1 — > e + e~ decay [30] at 
Tevatron, and the possibility to discover Z 1 gauge boson is analyzed in [31]. 

In the present work, we investigate the FCNC transition of — > in family 

non-universal Z' model. In particular, we analyze the differential branching ratio, forward- 
backward asymmetry as well as some related asymmetry parameters and polarizations 
and compare with the predictions of the SM. Note that the rare baryonic A;, — > A£ + £~ 
decay within family non-universal Z' model was analyzed in [32-34] within also family non- 
universal Z' model. The implications of non-universal Z' model on B meson decays were 
investigated in [35-43]. The B -> K 2 (-> K%)l+l~ [44], B -> K±l + l, and B -> K n [45,46] 
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were investigated in the same framework as well. The effects of a family non-universal Z' 
gauge boson is also searched for B — > im decays in [47]. 

The outline of the paper is as follows. In section 2, we present the effective Hamiltonian 
and transition matrix elements responsible for the E& — > transition. In section 3, 

we analyze the differential branching ratio, forward backward asymmetry, double lepton 
polarizations as well as some other related asymmetries in the Z' model and compare the 
obtained results with the SM predictions. The last section encompasses our concluding 
remarks. 

2 The — > transition in Family Non-universal 

Z' Model 



2.1 The Effective Hamiltonian 



Neglecting terms proportional to 



0(1(T 2 ), the effective Hamiltonian of the FCNC 



transition of E& — > H£ + £ , proceed via quark level b — > s£ + £ in the SM, can be written as 
[48-51] 



n 



eff 



C e 9 ff s lfl (l - 75 )&*y^ + C 10 s 7 m(1 - l^^ld 



- 2m h C e 7 ff ^sia liV q»{l+ lb )b£ 1 n 



q 



(2.1) 



where a em is the fine structure constant at Z mass scale, Gp is the Fermi coupling constant, 



Vij are elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and Cy , Cg^ , C 
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are the Wilson coefficients. When Z' boson is considered and Z — Z' mixing is neglected, the 
extra part which should be added to the above effective Hamiltonian is written as [35,36] 



H, 



eff 



2G F 
72 



n sb D U 

v tb v t * s 



S7 M (1 - 75)6*7^(1 - 7 5 )£ + 



BkBR 



u - 



V tb V t * s 



s 7m (1- 75 )6%(1 + 75K 



(2.2) 



where B^ b 



\B^ b \e^i and Bj 



L R correspond to the chiral Z' couplings to quarks and leptons, 
respectively. Hence, in Z' model, to get the effective Hamiltonian for the transition under 
consideration, we need to make the following replacements in Eq.(l) to consider Z' boson 
contributions besides the Z boson: 



Cf 



-> cf = cf 



n sb 1 D L 



47T 

— (28.82) 
a s K ; V tb V t * s 

C10 ->■ C[ = C10 H (28.82) sb [B\ t - B%), 

« s VtbVts 



(2.3) 



where a s is the strong coupling constant. Here we should mention that the Wilson coefficient 
C7 remains unchanged. 
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In SM, the Wilson coefficient C^* in leading logarithm approximation is written as (see 
[52]) 

C e / f ((i b ) = v ¥i C 7 (fi w ) + | (tP - 77 s ) C 8 (// w ) + C 2 (fi w ) ]T , (2.4) 

6 i=i 

where 

C 2 ( m ) = 1 , C 7 (/^) = -^o(^) , C a (fi w ) = ~E (x t ) . (2.5) 

The functions, D (x t ) and E (x t ) are defined as 

(8^ + 5^-7^) gg(2-3x t ) 

" 12(1 -x t ) 3 + 2(1 - Xt y lnxt ' (2 ' 6) 

^ = ii^ 2) + 2&^- (2 - 7) 

2 

where x t = -4?$- with m 4 and M w being the top quark and W boson masses, respectively. 

1V1 W 

The coefficients a, and hi are given as [53, 54] 



Oj = ( ||, ||, ^, -g, 0.4086, -0.4230, -0.8994, 0.1456 ), 

hi = ( 2.2996, -1.0880, -f , -0.6494, -0.0380, -0.0186, -0.0057 ). 



(2.8) 

3<±, — u.uoou, -u.uiou, -u.uuoi j. 

The parameter r\ in Eq. (2.4) is also defined as 



V = ^4 , (2-9) 



with 



«.(*) = > (2-10) 

l-/3 ^ln(^)' 
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where a s (m^) = 0.118 and /3 — , 

The Wilson coefficient Cg^(s') is written as [53,54]: 

C 9 e// (s') = Cg DR r}(s') + h(z, s') (3Ci + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) 
-^MM') (4C 3 + 4C 4 + 3C 5 + C 6 ) 

~h(0, s') (C 3 + 3C 4 ) + ^ (3C 3 + C 4 + 3C 5 + C 6 ) , (2.11) 
where s' — ^ with 4m^ < q 2 < (mj t — m s ) 2 , and 



^ = p ™ + ^)__ 4Zo(a . t) + p^ (a . t)j (212) 
sin i7 

3 



here, NDR stands for the naive dimensional regular izat ion scheme. We ignore the last term 
in this equation due to the negligible value of Pe- The Pq DR = 2.60 ± 0.25 [53, 54] and the 
functions Y (x t ) and Z (x t ) are defined in the following form: 



Yo(x t ) = | 



x t - 4 3x t 

7 + 7 7v7 In^t 

^ - 1 (x t - l) 2 



(2.13) 



and 



„ . . 18x 4 t - 163a; 3 + 259x 2 t - 108x t 
UXt) = U4(x t - 1)3 + 



32xf - 38xf - 15x 2 + 18x t 1 



72(x t - l) 4 



In Eq.(2.11), the r](s') is given as 

rj(s') 

with 



IT 



2s'(l + s')(l-2s'), -/ 5 + 9s'-6s' 2 
■ In s + 



3(l-s / ) 2 (l + 2s / ) 6(l-s / )(l + 2s')" 

At /it scale, for the coefficients Cj (j = 1, ...6) we have 



\nx t . 
(2.14) 

(2.15) 



(2.16) 



Cj — ^2 kjif] ai , 



(2.17) 



where kji are given as: 



hi = ( 


o, 


o, 


1 

2' 


1 

2' 


o, 


o, 


o, 


o ), 


= ( 


o, 


o, 


1 

2' 


1 

2' 


o, 


o, 


o, 


o ), 


hi = ( 


o, 


o, 


1 

14' 


1 

6' 


0.0510, 


-0.1403, 


-0.0113, 


0.0054 ), 


/C4j = 1 


o, 


o, 


1 

14' 


1 

6' 


0.0984, 


0.1214, 


0.0156, 


0.0026 ), 


hi = 1 


o, 


o, 


o, 


o, 


-0.0397, 


0.0117, 


-0.0025, 


0.0304 ), 


hi — 1 


o, 


o, 


o, 


o, 


0.0335, 


0.0239, 


-0.0462, 


-0.0112 ). 



(2.18) 



The other functions in Eq. (2.11) are also given as: 
h(y,s') = 



8, m b 8, 8 4 

■- In my H 1 — x 

9 n b 9 y 27 9 



--{2 + x)\l-x\^ 



In 



l-z-i 
2 arctan 1 



— m 



X-l' 



for x 
for x 



s' 

s' ' 



(2.19) 
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where y — lory — z — — and, 



We also consider the long distance contribution (Y LD ) coming from J ftp family resonances 
to Cg^ and parameterized using Breit-Wigner ansatz [48] as 

Yld - ^h J^^S 1 ^ ■ P-21) 



where C^ - 1 = 0.362, and we consider only two lowest resonances J/ip(lS) and ip(2S) and 
choose the corresponding phenomenological factors K\ — 1 and «2 = 2. We use the exper- 
imental results on the masses and total decay rates of dilepton decays of the considered 
vector charmonium states [55]. For more details about the calculation of long distance 
contributions, see for instance [56,57]. 

Finally, the explicit expression for C w is given as: 



C = -£&L , (2.22) 
sin t>w 



where, sin 2 9 W = 0.23. 



2.2 Transition Matrix Elements and Form Factors 

The £5 — > T,fi + fi~ decay amplitude is obtained by sandwiching the effective Hamiltonian 
between the initial and final states 

M = (E(p)|tf eff |E 6 (p + <?)>, (2.23) 

which reads 

M = ^^^ {cf (E(p)| s - 7/i (l- 75 )fe|E & (p + ,))^ 

+C[ (E(p)\s^(l-^)b\E b (p + q))h^ 5 l 
-2m b Cf i (E(p)|^^^(l + 75)&|S b (p + g))ZVz| . (2.24) 

To calculate the amplitude, we need to parameterize the transition matrix elements (S(p) |s7 M (l- 
75)6|Eb(p + g)) and {T l {p)\sa liv q v {l + 7 5 )6|E&(p + g)) in terms of twelve form factors fi, gi, 
ff and gj (i = 1, 2, 3) as follows: 

(E(p) I s 7m (1 - 75 )& I S 6 (p + q)) = u E (p) [iMv 2 ) + ™^ V h{q 2 ) + q^fsiq 2 ) 

- i^giiq 2 ) - iv»viM u g2{q 2 ) - q^i593(q 2 )]u^ b (p + q) , 

(E(p) I sia,„q»(l + 75 )& I E 6 (p + q)) = « E (p) ^(g 2 ) + i°^q v fl{q 2 ) + ^/ 3 V) 

+ 7m753iV) + i^l 5 q u 9l(q 2 ) + g'W \^)W b (p + ?) , ( 2 -25) 
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where u-£ b and «e are spinors of £5 and £ baryons, respectively. 

The form factors as the main inputs in the analysis of the £& — > s Lt r t~ have been very 
recently calculated in full theory via light cone QCD sum rules in [3]. The fit function of 
transition form factors is given as 

/ftfJbPV)] = t-^Vt + 7j — 3?Lw' < 2 ' 26 > 

where the fit parameters a, b, and mf it are presented in Table 1. 





a 


b 


rrifu 


q 2 = 


h 


-0.035 


0.13 


5.1 


0.095 ±0.017 


h 


0.026 


-0.081 


5.2 


-0.055 ±0.012 


h 


0.013 


-0.065 


5.3 


-0.052 ±0.016 


9i 


-0.031 


0.15 


5.3 


0.12 ±0.03 


92 


0.015 


-0.040 


5.3 


-0.025 ±0.008 


93 


0.012 


-0.047 


5.4 


-0.035 ±0.009 


ft 


1.0 


-1.0 


5.4 


0.0 ±0.0 


ft 


-0.29 


0.42 


5.4 


0.13 ±0.04 


ft 


-0.24 


0.41 


5.4 


0.17 ±0.05 


9l 


0.45 


-0.46 


5.4 


-0.010 ±0.003 


T 
92 


0.031 


0.055 


5.4 


0.086 ± 0.024 


T 
93 


-0.011 


-0.18 


5.4 


-0.19 ±0.06 



Table 1: Parameters appearing in the fit function of the form factors, f\, f 2 , ^3, gi, #2, 93, 
fT, ft, ft, 9i, 92 an d gt in full theory for £5 — > Y,l + l~ together with the values of the 
form factors at q 2 = [3]. 
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3 Observables Related to the £& — > YjI + 1 Transition 



3.1 Branching Ratio 

Using the amplitude mentioned above, the differential decay rate is found as (see also [58]) 



dF(z, s) _ G 2 F a 2 m m^ 
ds dz ' 16384tt 5 



■\v tb v:fvV\ 



%(s)+T 1 (s)z + T 2 (s)z 2 



(3.27) 



where s = and z = cos 9 with 9 being the angle between the momenta of E& and £ + in 
the center of mass of leptons, A = A(l, r,s) = 1 + r 2 + s 2 — 2r — 2s — 2r§, r = m\jm\ and 



v = x 1 - 



Am 2 . 



The functions, %(s), 7i(s) and T 2 (s) are given as 



%{s) = Z2m 2 t mi h s(l +r-s) (\D 3 \ 2 + \E 3 \ 2 ) 
+ 64mfm| 6 (l - r - s) Re[D{E 3 + D 3 E{] 
+ 64ml b Vr{Gm 2 e - m^Re^^i] 

+ 64m 2 m| 6 v / ^(2m Ei) sRe[ J D^3] + (1 - r + s)Re[£>^ 3 + ^s]) 
+ 32m| b (2mf + m| fc s){(l - r + s)m S6V ^Re[A* 1 A 2 + BJB 2 ] 
- m S6 (l - r - s) Re\A\B 2 + A^i] - 2^(Re\A\B l ] + m| 6 s Re[^S 2 
+ 8m| fc {4m 2 (l + r - s) + m\ h [(1 - r) 2 - s 2 ] } (|Ax| 2 + |#i| 2 



where 



+ 8m| 6 {4m 2 [s + (1 + r - s)s] + m| 6 s[(l - r) 2 - s 2 ] } (|A 2 | 2 + |5 2 | 2 ) 

- 8m| 6 {4m^(l + r - s) - m\ b [(1 - r) 2 - s 2 ] } (|Di| 2 + |£i| 2 ) 

+ 8m| b st; 2 { - 8m Sb s^Re[D* 2 E 2 ] + 4(1 - r + s)^Re[D\D 2 + E{E 2 ] 

- 4(1 - r - s) Re[£>If? 2 + L^] + m Efc [(1 - r) 2 - s 2 ] (|Z^ 2 | 2 + |£ 2 | 2 ) }, 



(3.28) 



%_{s) = -16m 4 At svV^{2Re(A* 1 D 1 ) - 2Re(B{E l ) 

+ 2m A6 i?e(5* J D 2 -5 2 * J D 1 + A; J B 1 -A*E 2 )} 

+ 32m 5 Ab s vVz{m Ab (l - r)Re(A* 2 D 2 - B* 2 E 2 ) 

+ y/rRe(A* 2 Di + A*/J 2 - B\E X - B\E 2 )) , 



T 2 {s) = -8m|yS (l^l 2 + Ifixl 2 + |A| 2 + |#i 

+ 8m|^ 2 £( |A 2 | 2 + |E 2 | 2 + \D 2 \ 2 + |£ 2 | 2 ) , 



(3.29) 



(3.30) 



A 1 = -'^r ( f-i- 



C e 7 ff (fl + gl) + CI" (/1 - , 
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A 2 


= A 1 (1^2), 


A 3 


= M (1 ->■ 3) , 


B t 


= Ai (ft -> ft T ->■ - 


D 1 


= Cio (/i - ft) , 


D 2 


= A(1^2) , 


Dz 


= £>i(1^3) , 


Ei 


= A ->■ • 



To numerically analyze the differential branching ratio with respect to q 2 , we perform 
integral over z in the interval z G [—1,1] in Eq. (3.27) and take the values of input 
parameters as m t = 167 GeV, mw = 80.4 GeV, mz = 91 GeV, m& = 4.8 GeV, m c = 
1.46 GeV, fi b = 5 GeV, fi w = 80.4 GeV, \V tb V*\ = 0.041, G F = 1.17 x 10~ 5 GeV" 2 , 
a em = ^, m s = 1.192 GeF, m Efc = 5.807 GeV and m M = 0.1056 GeV [55]. The other 
input parameters related to the family non-universal Z' model are given in Table 2, where 
SI and S2 are two sets of parameters correspond to UT-fit Collaboration results [59]. These 
two different sets are obtained fitting the results to B s — B s mixing data [59]. The results 
of the works [45,46] indicate that 52 is more favored than SI. 

In order to maximize the effect of Z' boson, we consider the maximum values of the 
related parameters, i.e. \B%,\ = 1.31 x 10~ 3 , tp^ = -79°, Bft + Bft = -6.7 x 10" 2 , and 
Bji - B* = -9.3 x 10- 2 for SI set and \B%\ = 2.35 x lO" 3 , tf® = -86°, B^ + fl* = 
-2.60 x 10" 2 , and - B* = -3.4 x 10" 2 for S2 case. 









B i» + B w 


f>L tdR 

MM MM 


SI 


(1.09 ±0.22) x 10- 3 


-72 ±7 


(-2.80 ±3.90) x 10~ 2 


(-6.70 ±2.60) x 10~ 2 


S2 


(2.20 ±0.15) x 10~ 3 


"82 ±4 


(-1.20 ± 1.40) x 10" 2 


(-2.50 ±0.90) x 10~ 2 



Table 2: The values of the input parameters for the Z' couplings in two sets. 

We plot the differential branching ratio, forward backward asymmetry and other asym- 
metry parameters with respect to q 2 to check the sensitivity of these observables to SI and 
S2 sets of parameters. For this aim, we plot the variation of the differential branching ratio 
in terms of q 2 in figure 1. From this figure, we obtain the following results: 

• we see considerable discrepancies between the Z 1 model and the SM predictions for 
SI set of parameters. In the case of S2, the difference is observed at higher values of 

• The order of differential branching ratio indicates the possibility of detecting this 
decay mode at LHC. 
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Figure 1: The q 2 dependence of the differential branching ratio for the S b — > decay 
for SI, and ,52 sets of parameters as well as the SM. 

3.2 Forward-backward Asymmetry 

One of the most promising tools in detecting the NP effects is the forward-backward asym- 
metry of leptons which is defined as 

N f -N b , . 

-4- = ^Ti- (3-32) 

where Nf is the number of events that particle is moving "forward" with respect to any 
chosen direction, and iVj, is the number of events that particle moves to "backward" direc- 
tion. In technique language, the forward-backward asymmetry Afb(s) is defined in terms 
of the differential decay rate as 

a /-\ Jo dsdz J-i dsdz / Q qq\ 

Afb(s) - dT , (3.33) 

Jo dMz^ s)dZ + LdMz iz - s)dZ 

We depict the dependence of Afb on q 2 in figure 2. From this figure, we conclude that 
the predictions of Z' model for SI and .52 sets deviate considerably from the SM result 
especially at higher values of q 2 . 



3.3 Baryon polarization asymmetry parameter - 

Asymmetry parameters characterize the angular dependence of differential decay width for 
the cascade decay S 6 — > £(— >■ a + b)V*(— > with polarized and unpolarized heavy 

baryons. In this part, sensitivity of the baryon polarization asymmetry parameter to new 
Wilson coefficients is analyzed. The helicity amplitudes for the — > ~E£ + £~ decay are 
obtained by analyzing quasi two body decay £5 — > T.V*, followed by the leptonic decay of 
V* — > in [60]. These amplitudes are obtained using helicity amplitude formalism and 
polarization density matrix method, demonstrated in [61,62]. Considering 

S fe 1/2+ -)• Y}' 2+ (^a + b) + V*(^ £ + T) , (3.34) 
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Figure 2: The q 2 dependence of the forward backward asymmetry for the S b — > 
decay for 51, and S2 sets of parameters as well as the SM. 



with V* being off-shell 7 or Z boson, the normalized joint angular decay distribution for 
the two cascade decay is written as [63-68] 



dT{q 2 ) 



dq cfcos 9 dcos 9% 



GpOL e 



■v tb v t 



Is 



J\(ml b ,ml,q 2 )^/\(ml,m 2 a ,m 2 b ) ,m.„2 

vB(h b — > a + 6) \M\ , 



10247r 3 m| m| 



(3.35) 



where M. is calculated in [60]. For the leptonic part, in the rest frame of the intermediate 
boson, the angle of the anti-lepton with respect to its helicity axes is shown by 9. For the 
hadronic decay, in the rest frame of the £&, 9^, is the angle of the a momentum with respect 
to its helicity axes. 

The polar angle distribution of £ — > a + b decay is obtained by integrating Eq. (3.35) 
with respect to 9 and it gives the differential decay rate in terms of q 2 and 9% [60] 



dT{q 2 ) 



1 + a(q 2 )a^(q 2 ) cos6» s (g 2 ) , 



dq 2 d cos 9y, 

where the baryon polarization asymmetry parameter a^, is given as 



(3.36) 



«s(g 2 ) 



-{4m 



A 



+1/2.+1 



3A(q 2 
- 2m 2 |A+i/2,o| -g 2 ( 

-4-1/2 



+ 2g 2 ( 



1+1/2,+! 



l +l/2,0 



4m, 2 



2 



5 



+ w 
2 



B 



+1A+1 



+1/2,0 



) - 6mj; 



+l/2,t 



-1/2,-1 



1/2,-1 









2 2 










-4-1/2,0 




-4-1/2,0 


+ W 2 


£-1/2,0 


2 ) +6m 2 


B-l/2,t 


') 



+ 2m, 

and the asymmetry parameter a is obtained as 

2 



(3.37) 



a(g 2 ) 



3A(g 2 



-|4m 



l +l/2,+l 



+ 2q 2 ( 
10 



H-1/2,+1 



B 



+1/2.+1 



+ 2mj 
- Ami 

~q 2 



l +l/2,0 



+ 6m, 



B 



+l/2,t 



+ <?{ 



+ w 



"4+1/2,0 
2 



B-l/2, 



+ v 

2 

-1 



5 



+1/2,0 



1/2,0 



5_ 



1/2,0 



■ \ 2 

J - 2m^ A_i/ 2)0 



— 6m„ 



B-i/2,t\ 2 } , (3.38) 



with 



A(g 2 ) = |{4m, 2 



A. 



+1/2.+1 



+ 2g 2 ( 



A 



+ 2m] U+i/2,1 



+ 6m/, 



5 



+1/2.+1 
2 



5 



+l/2,t 



+ <? 2 ( -4+1/2,0 



+1/2.+1 
2 



4mi 



"4-1/2,-1 

. 2 
-4-1/2,0 



V( 1-4-1/2, 



-B-l/2, - 



+ V 

2 

1 



5 



+1/2,0 



B_ 



1/2,0 



) + 2m 2 



1/2,0 



+ 6m» 



i/2,*| 2 }. (3.39) 



The definitions for A\^\ v and B\ u \ v are given in [60] where Aj and \y are helicities of lepton 
pairs and vector boson, respectively. We plot the dependence of the baryon asymmetry 
parameter as on q 2 in figure 3. From this figure, it is clear that 
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Figure 3: The q 2 dependence of the asymmetry parameter a s . 



• in general, the results obtained using two sets of parameters deviate significantly from 
the SM predictions. 

• The effect of Z' boson on baryon asymmetry parameter for S b — > T,fi + fi~ channel is 
larger compared to that of A b — > A/i + /i~ decay mode analyzed in [33]. 



3.4 Polar angle distribution parameters- cte and f3e 

The polar angle distribution of V* — > l + l~ decay is obtained by integrating Eq. (3.35) with 
respect to #s which gives the differential decay rate in terms of q 2 and 9 



dT{q 2 



dq 2 d cos 9 



~ 1 + 2a e cos 9 + (3g cos 2 9 , 



(3.40) 
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where 



and 



a e (q 2 ) = 



Ai(g 2 



2vq 2 Re A +1/2 ,+iB* +1/2t+1 - A_ 1/2 ,-1-8*1/2,-1 



(3.41) 



Ax 

+ 4m 2 

- 4m 2 

+ Am 2 p 



A 



+1/2, +1 



+ g 2 ( 



A 



+i/2,+i 



B 



+1/2.+1 



l +l/2,0 



-1/2,-1 



1/2,0 



2 + 4 



l +l/2,0 



2 



5 



+1/2,0 



) 



-1/2,-1 
2 

/2,0 



+ V 



+ V 
2 



-1/2,-1 
2 



B 

#-l/2,o| )} , 



) 



(3.42) 



with 



Ai(g 2 ) = 4m 2 



+ 4m 2 ; 
+ 4m 2 
+ 4m» 



5 



+ l/2,+l 
+1/2,*' 



^-1/2, 



l/2,t 



2 + « 2 ( 
+ 9 2 ( 

2 2 

+ g 2 ( 



A 



+1/2.+1 
2 



2 \ 



l +l/2,0 



2 



5 



+1/2,0 



1/2,-1 
2 

-4-1/2,0 



5_ 



2 \ 

#-1/2,-1 J 



(3.43) 



We plot the dependence of the polar angle distribution parameters ag and j3g on q in figures 
4 and 5, respectively. From these figures, we read that 




Figure 4: The q 2 dependence of the asymmetry parameter ag. 



• similar to the A& — > A/i+yU - channel considered in [33], the asymmetry parameter f3g 
in Z' model gives roughly the same result as the SM; however in the case of ag there 
is a considerable discrepancy between the two model predictions. 

• The zero point of ag parameter is shifted to the right for both sets of parameters 
compared to the SM prediction. 

• The effect of Z' boson on the ag parameter in E b — > channel is large compared 
to the A b — > Afi + fi~ channel discussed in [33] . 
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Figure 5: The q 2 



dependence of the asymmetry parameter (3g. 



3.5 Double Lepton Polarizations 

To analyze the double lepton polarizations, we need to define the following orthogonal unit 
vectors sf' 1 with % = L, T or N in the rest frame of double leptons (see for instance [69-71] 
for detailed information) 



where p± and p-^ are the three-momenta of the leptons £ ± and S baryon. The symbols L, 
T and N denote longitudinal, transverse, and normal polarizations, respectively. Using the 
Lorentz boost, these unit vectors are transformed from the rest frame of the leptons into 
the center of mass (CM) frame of them along the longitudinal direction. As a result, we 
get 



where p + = —p_ and En and mi are the energy and mass of leptons in the CM frame, 
respectively. Under the above transformation, the other two unit vectors, s^, remain 




(3.44) 




(3.45) 
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unchanged. The double lepton polarizations are given as [69-71] 



/ dT(Sj ,s+) _ dTj-Sj ,s/) n _ / drjsj , -s+) _ dTj-Sj , -s/) \ 

V d,S d,S ' v <7.«. <7.s / 



/ dT{Sj ,Sj ) gX-g ,Sj k / ui ^ , -a, ) m y-Sj , —ay ) \ 



(3.46) 



The first (second) subindex of represents polarization of lepton (anti-lepton). Using the 
above definitions, some double lepton polarizations are obtained as 

167rmi m^vA f . 

PLiv(S) = -PiVL(s) = ^ * Im | (l _ + 

+ m Sb s(A*E 3 - + 5*L> 3 - 

+ mz b V?s(AlD 3 + + 5*£ 3 + B* 2 E X ) - ml b s 2 (B* 2 E 3 + A* 2 D 3 )^ , (3.47) 

PMs) = 16 ^|y Re{(l - r) ( l^l 2 + i^l 2 ) - s(A lDl - Bl E t ) 



B 1 D* + (A 2 + D 2 - D 3 )E{ - A ± E* - (B 2 -E 2 + E 3 )D\ 



+ ml b s(l-f)(A 2 D*-B 2 E*) 



+ m^Vf-s^D; + (A 2 + D 2 + D 3 )D\ - B X E* 2 - (B 2 - E 2 - E 3 )El 
-ml b s 2 (D 2 D* 3 + E 2 E*)} , 



(3.48) 



P NT (s) = -Ptn = 3A / ( g) lm \( A i D *i + B i E *i) + ml b s(A* 2 D 2 + B* 2 E 2 ) \ , (3.49) 



NN 



'§) = 



32m f 

3sA'5) Re | 24 ^t^K A i B i + D i E i) - 12m S6 m^v / ?s(l - f + s)(Ai^ + fii^) 



+ 6m S6 m|s[m Sb s(l + f - s)( |D 3 | 2 + |£ 3 | 2 ) + 2v^(l - r + sOPi^s + E x El) 
+ 12m Si m 2 s(l - f - s)(A lB * 2 + A 2 £* + + A^D 

- [As + 2mf(l + f 2 - 2f + fs + s - 2s 2 )]( |Ax| 2 + l^] 2 - |Dx| 2 - |£x| 2 ) 
+ 24m| ft m 2 v / fs 2 (A 2 5 2 * + D 3 E* 3 ) - ml b Xs 2 v 2 (\D 2 \ 2 + \E 2 \ 2 ) 

+ m| 6 s{As-2mf[2(l+r 2 ) - s(l + s) - r{4 + s)]}(\A 2 \ 2 + |5 2 | 2 )| , (3.50) 



32m 4 f 

Ptt(s) = ^§)H - 24m 2 v / fs(A lJ B 1 * + D 1 E^) - l2m Sb m 2 V?s(l - f + s)(£>i^ + E^) 
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- 24ml b m 2 e Vfs 2 (A 2 B* 2 + D 3 E* 3 ) 

- 6mx b m 2 e s[mj; b s(l + r - s)(\D 3 \ 2 + \E 3 \ 2 ) - 2<s/r(l ~ ? + s)(A ± A* 2 + B ± B*) 

- 12m s > 2 s(l - f - s)(A-iBl + A 2 B\ + D ± E* 3 + D 3 E*) 

- [\s-2ml(l + f 2 -2f + fs + s-2s 2 )]( y \A 1 \ 2 + \B l \ 2 ) 

+ mlj{Xs + m 2 [4(l - ff - 25(1 + r) - 2s 2 ]} ( \A 2 \ 2 + \B 2 \ 2 ) 
+ {As - 2m 2 [5(l - f) 2 - 7s(l + f) + 2s 2 ]}( l^l 2 + l^ 2 

- ml\s 2 v 2 [\D 2 \ 2 + \E 2 \ 



(3.51) 



where 



a'(s) = To(s) + ir 2 (s). 



(3.52) 



We plot the dependence of the double lepton polarizations on q 2 in figures 6 — 10. 
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Figure 6: The q 2 dependence of Longitudinal-normal (LN) polarization of leptons 




Figure 7: The same as figure 6, but for LT polarization 
From these figures, we observe that 
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Figure 8: The same as figure 6, but for TN polarization 





Figure 10: The same as figure 6, but for TT polarization 



• in general, there are considerable discrepancies between the Z' model and the SM 
predictions on the double lepton polarization asymmetries under consideration for 
both sets of parameters except for the Plt- 

• The big discrepancies of the Z' model results from the SM predictions appear at lower 
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values of q 2 for both sets. 



• In the case of Pln and Ptn, the effect of Z' boson is large for 52 set of parameters. 
For the P NN and P T t, the SI set of parameters is more effective than the S2 at lower 
values of q 2 . 

Determination of zero points of double lepton polarization asymmetries and comparison of 
the obtained results with the theoretical predictions can give valuable information about 
the existence of the Z' boson. 

4 Conclusion 

In the present study, we analyzed the E& — > EyU + /i~ decay channel in both the SM and 
Z' model with two sets of parameters. Our results on the considered physical quantities 
overall depict considerable discrepancies between the Z' and SM models predictions. These 
discrepancies can be considered as signals for existing the extra Z' gauge boson. The 
order of the differential branching ratio indicates that this channel can be studied at LHC 
in near future. Any measurement on the physical quantities considered in the present 
work and comparison of the obtained data with the theoretical predictions can give useful 
information not only about the existence of the Z' gauge boson but also about the nature 
of participating baryons, E 6 and E. Such comparison will also help us to understand which 
set of parameters is favored. Our numerical calculations lead that for some observables 
like forward-backward asymmetry, baryon polarization asymmetry parameter, polar angle 
distribution parameters; and Pnn and Ptt from double lepton polarization asymmetries, 
the SI set is more effective. While for differential branching ratio; and Pln and Ptn as 
other double lepton polarizations, the S2 set is favored. Similar analysis has been made for 
Aft — )> Afi + fi~ decay channel in [33] which is in agenda of different experiments nowadays. 
The study of different FCNC channels can provide us with more data which may help us 
in searching for Z 1 gauge boson as new physics effect. 
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